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Abstract: We present the derivation of the NNLO two-parton final state contributions 
to top pair production in the quark-antiquark channel proportionnal to the leading colour 
factor N^. Together with the three and four-parton NNLO contributions presented in a 
previous publication, this enables us to complete the phenomenologically most important 
NNLO corrections to top pair hadro-production in this channel. We derive this two-parton 
contribution using the massive extension of the NNLO antenna subtraction formalism and 
implement those corrections in a parton-level event generator providing full kinematical 
information on all final state particles. In addition, we also derive the heavy quark con¬ 
tributions proportional to N^- Combining the new leading-colour and heavy quark contri¬ 
butions together with the light quark contributions derived previously, we present NNLO 
differential distributions for LHC and Tevatron. We also compute the differential top quark 
forward-backward asymmetry at Tevatron and find that our results are in good agreement 
with the measurements by the DO collaboration. 
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1 Introduction 


The precise study of top quark properties provides a detailed probe of the Standard Model 
of particle physics and models beyond it. Due to its small lifetime, the top quark decays 
before it hadronizes and it is the only quark whose production dynamics can be studied 
without having to account for hadronization effects. Following its discovery twenty years 
ago [1, 2], the top quark has been studied intensively at the Tevatron. Owing to the 
limited number of top quarks produced, differential studies there suffered from rather large 
uncertainties [3, 4]. The measured top quark forward-backward asymmetry received much 
attention since it was first found to deviate substantially from Standard Model expectations 
[5]. 

With the large number of top quark pairs produced at the LHC, the study of its prop¬ 
erties has become precision physics, and the measurement of very precise differential cross 
sections, attainable. Recently, the ATLAS and CMS collaborations reported measurements 
of normalised differential observables in tt production in several kinematical variables, snch 
as the transverse momentum, rapidity and invariant mass of the tt system, as well as the top 
quark transverse momentum and rapidity [6-10]. More recently, both LHC collaborations 
have measured the charge asymmetry in top-pair production, [11, 12]. 

These measurements will allow for a very detailed and accurate probe of the top quark 
production mechanism. To reliably interpret data, these very precise measurements must be 
matched onto equally accurate theoretical predictions, which can be obtained by including 
contributions up to next-to-next-to leading order (NNLO) in perturbative QCD. At present, 
NNLO corrections have been included in calcnlations of the inclusive total tt cross section 
in [13], for the inclusive and differential Tevatron top quark forward-backward asymmetry 
in [14], and for differential LHC distributions in [15]. 

At NNLO, perturbative calculations of collider observables are typically carried out 
using parton-level event generators. These programs generate events for all parton-level 
subprocesses relevant to a given final state configuration up to NNLO accuracy and pro¬ 
vide full kinematical information on an event-by-event basis. An NNLO event generator 
for observables with n final state particles or jets involves three main building blocks: the 
two-loop corrections to the n-parton final state, referred to as double-virtual contributions 
(dcr^Jj^LQ), the one-loop corrections to the (n-|- l)-parton final state, called real-virtual con¬ 
tributions (du^NLo)’ tree-level {n + 2)-parton double real contribution (du^^LQ). 

These three building blocks involve infrared divergences that arise from the exchange or 
emission of soft and collinear partons and cancel only in their sum. In addition, two mass 
factorisation counter-terms, du^j^Lo and dcr^j^^Q, are needed in the three and two-parton 
final states contributions respectively in order to cancel infrared divergences originated 
from initial state collinear radiation. 

The combination of snbprocesses of different particle multiplicity which are individually 
infrared divergent is a major challenge in the construction of NNLO parton-level event 
generators. Employing a subtraction method to regulate these infrared singnlarities, the 
NNLO partonic cross section for top pair production in a given partonic channel has the 


- 2 - 


general structure [16] 
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Two types of subtraction terms are introduced: dd^j.^Lo 4-parton final state, and 

ddJJ^LO 3-parton final state. The former approximates the infrared behaviour of the 

double real contributions dd^^^Q in their single and double unresolved limits, whereas the 
latter reproduces the single unresolved behaviour of the mixed real-virtual contributions 

dd^v 

‘^'^NNLO- 

In the context of the present computation, we further decompose the double real 
subtraction term dd^j.^LQ as in [17], which allows us to rearrange the different terms in 
eq.(l.l) into the more convenient form 
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The NNLO contributions to the top pair production in the quark-antiquark channel 
can be decomposed into colour factors as 
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with Nc being the number of colours, Ni the number of light quark flavours, and the 
number of heavy quark flavours. All coefficients multiplying the different colour factors in 
eq.(1.5) are gauge invariant and can be computed independently. 

The goals of this paper are twofold. The first aim is the derivation of the two-parton 
contribution given in eq.(1.2) for the colour factor within the theoretical framework of 
antenna subtraction with massive fermions. This framework has been established in the 
massless case in [16, 18-23] and extended to the massive case in [17, 24-29]. We shall not 
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repeat it here. Together with the three and four-parton contributions derived in [17] the 
two-parton contribution yields the coefficient A in eq.(1.5) and enables us to obtain theoret¬ 
ical predictions for the most significant NNLO contribution to top pair hadro-production in 
the quark-antiquark channel. In addition, we derive the heavy flavour contributions labeled 
in eq.(1.5) as and G^, presenting a detailed discussion of the ultraviolet renormalisation 
of the loop amplitudes involved. 

Our second goal is to present phenomenological results for Tevatron and LHC energies, 
including the leading-colour and heavy quark contributions derived in this paper, as well 
as the light quark terms, which we completed in previous papers [28, 29]. 

We employ analytic expressions for all matrix elements, subtraction terms and inte¬ 
grated subtraction terms that arise in our calculation. For the real-virtual contributions, 
we use DpenLoops [30] in combination with CutTools [31] and for the double virtual con¬ 
tributions we use the analytic two-loop matrix-elements of [32, 33].^ As a result, at the 
real-virtual and virtual-virtual levels, the explicit infrared poles are cancelled analytically. 
Furthermore, the numerical evaluation of the finite remainders at each level included in 
the parton-level generator is substantially faster than when these terms are evaluated nu¬ 
merically. 

The paper is organized as follows: In section 2 we present the derivation of the new in¬ 
tegrated massive initial-final antennae required in the double virtual counter-term d(Tjjj.,fLo 
for the leading-colour contributions to top pair production in the quark-antiquark channel. 
Section 3 contains the explicit derivation of the counter-term ddJjj.^LQ expressed in terms of 
massive integrated dipoles. Employing those dipoles, we find that the virtual-virtual sub¬ 
traction term has a similar structure to the one found in the massless case in the context 
of di-jet production at NNLO [18]. Section 4 is dedicated to the heavy quark contribu¬ 
tions proportional to and Nh/Nc- Section 5 contains our numerical results: Various 

differential distributions are presented and the phenomenological impact of this computa¬ 
tion is discussed. Finally section 6 contains our conclusions. We enclose three appendices 
containing the master integrals required to compute the integrated antennae presented in 
section 2, as well as the phase space parametrisations needed to calculate these integrals. 


2 Integrated initial-final massive antennae 
2.1 General features 


One of the aims of this paper is the computation of the NNLO two-parton contribution 
proportional to for top pair production in the quark-antiquark channel. Within the 
antenna subtraction method, this involves the construction of the double virtual counter¬ 
term da^- which renders the two-parton contribution 
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finite and suitable to be implemented in an NNLO parton-level generator. 

^For the gluon-gluon channel, which is beyond the scope of the present paper, analytic expressions for 
the leading-colour and light quark contributions of the two-loop amplitudes have been presented in [34, 35]. 
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As shown in eq.(1.4), this virtual-virtual counter-term contains mass factorisation 
terms together with integrated double real and real-virtual subtraction terms. In order 
to compute these, the integrated forms of the antenna functions present in the correspond¬ 
ing real-virtual and double real subtraction terms are needed. The integrated antennae 
are obtained by integrating the antenna functions over the corresponding antenna phase 
space inclusively. This integration is carried out analytically, and turns the implicit soft 
and collinear singularities of the antenna functions into explicit poles in the dimensional 
regularisation parameter e, which are then cancelled against the explicit poles of the double 
virtual matrix elements and mass factorisation counter-terms. 

In the context of this paper, we need the integrated forms of two massive initial-final 
A-type antennae which are presented here for the first time. These antennae involve a 
massless initial state quark and a massive quark in the final state as radiators, and one 
or two final state gluons which are unresolved, i.e. soft and/or collinear. More concretely, 
the calculation that we present in this paper requires the integrated forms of the tree- 
level four-parton antenna 3^, 4^, 2g) and the one-loop leading-colour three-parton 

antenna Ag’'^(I q, 3g, 2q). The tree-level four-parton antenna appears in the double real 
subtraction term (dd t ^ 2 ) given in [17]. It is needed in order to capture the infrared 
behaviour of the leading-colour double real contributions associated to the partonic process 
qq —7- tigg when the two final state gluons are unresolved and colour-connected. The one- 
loop antenna appears in the real-virtual subtraction term given in [17]. It is 

required to capture the infrared behaviour of the leading-colour part of one-loop matrix 
element associated to the partonic process qq —)• tig. 

The general definition of the integrated forms of three and four-parton massive anten¬ 
nae in all three configurations (final-final, initial-final, initial-initial), as well as the cor¬ 
responding phase space factorisations and mappings, have been presented in [24, 27, 28]. 
Below we will only briefly recall the definition of initial-final integrated massive antennae, 
as they are the only ones needed in the calculation presented in this paper. 

Initial-final three-parton tree-level and one-loop antennae like 3 c,, 2q) {I = 0,1) 

are generically denoted as with parton i in the initial state, and j,k in the final 

state. The kinematics associated to these initial-final antenna functions with one massive 
radiator are Pi + q ^ Pj + Pk, with q^ = —Q^ < 0, = 0 and = tuq, and the 

integrated antennae are defined as [28] 



( 2 . 2 ) 


d<h 2 is the corresponding 2 —)• 2 phase space, and 



(2.3) 


The initial-final massive antenna phase space denoted by d<hx, is given by 



(2.4) 
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In addition to the phase space integration, one-loop antennae must be also integrated over 
the loop momentum. 

Initial-final four-parton antennae like A^{1q, 3^,4^, 2g) are generically denoted as 
with parton i in the initial state and j, k, I in the final state. The associated kinematics is 
Pi + Q ^ Pj + Pk + Ph with = —Q^, Pi = p] = Pk — ^ pf = wlq. The integrated 
forms of these four-parton initial-final antennae are obtained as 

1 ~t~ TfiP ^) r 

- 2-k ^ J (2.5) 


and the corresponding antenna phase space is given by 


-|- mg) 

^*^Xijki{Pj,Pk,pi]Pi,q) = - ^ - d<^3{pj,Pk,pr,Pi,q)- 


( 2 . 6 ) 


After the inclusive phase space integration, the three and four-parton integrated an¬ 
tennae are functions oi Q^, pi ■ q and tuq. We change the dependence on the latter two 


variables for 
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We also change the dependence on for and in three and four-parton inte¬ 

grated antennae respectively, where the bars and tildes are related to the initial-final phase 
space mappings that can be found, for example, in appendix B of [28]. The invariants are 
defined as Sij = 2pi ■ pj with the crossings explicitly performed on the momenta pi, pj, i.e. 
Pi > 0- In integrated subtraction terms, the mapped final state momenta pjk and pjki 
are relabeled in accordance to the labeling of momenta in the lower multiplicity final state, 
in such a way that the integrated antennae end up depending on invariants . This is the 
notation that we shall follow throughout this paper. 

In all cases, the inclusive phase space integration is carried out following the standard 
technique of reduction to master integrals using integration-by-parts identities (IBP) [36, 
37] . The computation of the master integrals is then carried out either directly to all-orders 
or iteratively, order by order in e, using differential equations techniques [38]. 

All the master integrals found contain multiplicative factors of the form (1 — Xi)~^^ 
which regulate soft endpoint singularities and should be kept unexpanded. All other terms 
can be safely expanded in e. As we explain in detail in appendix C, in some of the 
master integrals found in the reduction of the one-loop antenna of interest here, namely 
A3’^'^(1 q, 3g,2g), more than one e-power of (1 — Xi) is encountered. Therefore, the master 
integrals which we denote collectively by Iaixi,xo,e) can be most generally expressed as: 


Ia{xi,xo,e) = ^(1 - Xi)”" ''"i?),”)(xi,xo,e). 
n^m 


(2.9) 


The integers m and n are specific to each master integral; the functions i?i"'^(xj, xq, e) are 
regular as x^ —)■ 1 and can be calculated as Laurent series in e. 
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The integrated antennae collectively denoted by X(xi,xo,e) are linear combinations 
of master integrals with coefficients containing poles in e as well as in (1 — Xj). After the 
masters have been inserted into the integrated antennae, these take the form 


X{xi,xo,e) = {I - Xi) ^ '^^TZx{xi,XQ,e) 


( 2 . 10 ) 


where TZx{xi^XQ, e) is a regular function as x, —)• 1. The Laurent expansion of the singular 
factor (1 — is done in the form of distributions: 



( 2 . 11 ) 


m=0 


with 



( 2 . 12 ) 


It is worth noting that in the functions R^\xi,xo,e) and Tlx{xi,xo,e), which are 


regular as Xj —>■ 1, the massless limit xq —?> 1 cannot in general be safely taken. This is 
due to the presence of terms of the form log^(l — xq) = log^{mQ/{Q‘^ + rng)), which are 
expected and correspond to quasi-collinear limits of the antenna functions [24, 27, 39]. 

In order to perform the Laurent expansion of the integrated antennae, we distinguish 
two regions: a hard region where x^ 7 ^ 1 and a soft region where x, = 1. The highest order 
in e needed in the expansion of each master integral is determined by the e and x^-dependent 
coefficient that multiplies the integral in the integrated antenna. In the soft region, since 
the expansion in distributions generates an additional 1 /e factor, the functions are 
required to one order higher in e than in the hard region. Further details concerning the 
integration of A\{1q^ 3g, 4^, 2^) and ^^’^'^(Iq, 3^, 2^) will be presented below. 

2.2 The integrated tree-level antennae -^qQg and ~^gQgg 

The computation of dd^ nnlo n'^ requires the knowledge of the integrated tree-level three 
and four-parton massive initial-final antennae •^qQg and ~^qQgg- The former has been 
derived in [24], and only its pole part will be recalled below. The integrated flavour- 
violating four-parton antenna Qgg is new and will be presented here for the first time. 

2.2.1 The integrated three-parton antenna AgQg 
The pole part of the integrated antenna Ag Qg is given by 


Voles (A°Qg(e,Sij,Xi)) = -2lJj^^(e, sj^)(I(l - Xj) -h rjg^(xi), (2.13) 


where is a colour-order infrared singularity operator which has been presented in [24] 
and reads 


e^E 
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2r(l-e) 


(2.14) 








(2.15) 


and rgg^(x) is a colour-ordered splitting kernel given by 

rjg (a;) = (^^<5(1 - a:) +Po(x) ■ 

As can be seen from the equations above, -4° has its deepest pole at 1/e^. In the 
virtual-virtual subtraction term that will be presented in section 3, we will need products 
and convolutions of two of these integrated antennae, and we will therefore require its 
Laurent expansion up to order e^. The resulting expression is too long to be presented 
here. It is included in an ancillary Mathematica file attached to the arXiv submission of 
this paper. 

The integrated initial-final antennae depend explicitly on one momentum fraction Xi 
carried by the initial state parton. In order to combine these integrated forms with other 
integrated subtraction terms depending on both momentum fractions x\ and X 2 carried by 
the incoming quark-antiquark pair, it is useful to make all antennae explicitly depend on 
both xi and X 2 - We achieve this by introducing delta functions in the following way: 

~^q,Qgi(^,Slj,Xi,X2) = j dX 2 .A° Qg(e, Sg- , Xi)5(l - X 2 ) (2.16) 

■^q,Q9i<^,S2j,X2,Xl) = j dxlAl^Qg{e, , X2)5{1 - Xl) (2.17) 

The same prescription will be applied to Qgg and Ag’qg. 

2.2.2 The integrated four-parton antenna AgQgg 

To compute the integrated initial-final massive antenna A^g Qgg we start by expressing 
all phase space integrals as cuts of two-loop four-point functions with two off-shell legs in 
forward scattering kinematics [40] , and reduce these four-point functions to master integrals 
with LiteRed [41]. 

As mentioned above, we consider the following DIS-like kinematics 

P 2 + q^Pi+P3+P4 (2.18) 

with P 2 = Ps = Pi = 0, pI = mg, = —Q^ < 0. The denominators that we find are: 

Di = pI- ml 
D 2 = pI 
D3 = pI 
D 4 = {pi - qf 

D 5 = iP3+P‘if 
Dq = {P3-P2f 
Dj = {P4-P2f 

Dfi = {pi + P3f - ml (2.19) 

where Hi, D 2 and H 3 are the cut propagators. In the reduction procedure we impose 
momentum conservation and mass-shell conditions for the external legs, and we discard all 
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(e) /[4,7,8] 



(c) -^[4,8] 


Figure 1: Master integrals for Integrals are labelled according to the denominators 

involved (see eq.(2.19)). Bold (thin) lines are massive (massless). The double line in the 
external states represents the off-shell momentum q with = —Q^. The cut propagators 
are the ones intersected by the dashed lines. 


integrals that do not contain the propagators Di, D 2 or 1)3 with power -1. We find the six 
master integrals depicted in fig.l. Their explicit expressions, together with a parametrisa- 
tion of the phase space measure are presented in appendices A and B respectively. 

The integrals /[q], /[8 ] and /[ 4 ] are known [24, 25]. We computed the remaining master 
integrals, namely, ^[4,8]) -^[4,5,8] .f[4,7,8]j using differential equation techniques, and we 

present them here for the first time. Although the differential equations for these three 
integrals are coupled, they can be decoupled order by order in e, which enabled us to solve 
them iteratively. We fixed two of the three sets of undetermined integration constants by 
demanding that the integrals be regular as x —)• 1 after factoring out an appropriate power 
of (1 — x). In order to fix the remaining set of integration constants we calculated via direct 
evaluation the soft limit of /[ 4 ^ 7 ^ 8 ]) employed it as a boundary condition. The details 
of the computation of this soft limit can be found in appendix B. 

The highest order in e with which each master integral contributes up to the finite 
part of Qgg is specified in table 1. As mentioned before, in the soft region (x = 1), they 
must be expanded to one order higher than in the hard region (x 7 ^ 1 ). 

After expanding in distributions we find that the integrated four-parton massive an¬ 
tenna -^qQgg has a deepest pole at C>(l/e^). This is expected since at the unintegrated 
level A^{lQ,3g,4:g,2g) contains double soft and triple collinear singular limits, which are 
employed to reproduce the singular behaviour of the double real contributions associated 
to the partonic process qq —)• ttgg. The integrated antenna can be written in terms of 
harmonic polylogarithms (HPLs) with arguments Xj or xq and generalised harmonic poly¬ 
logarithms (GPLs) of argument xq and weights involving 1/xj. GPLs and HPLs appear 
with up to trascendentality three and four respectively. Given its length, the full expression 
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Integral 
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eO 
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-^[4,5,8] 
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e-3 

-^[4,7,81 
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Table 1: Deepest poles and highest order in e needed for each master integral in Qgg 


of this antenna up to 0{e^) is given in the ancillary file attached to the arXiv submission 
of this paper. The deepest poles read 
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( 2 . 20 ) 


2.3 The integrated one-loop three-parton antenna A}q^Qg 

The leading-colour one-loop antenna employed in [17] for the construction of a real-virtual 
subtraction term for partonic process qq —)• ttg at one-loop is given by 


^ 3 ’ (lQ)3g,2q) 


d°(l 3 2 ) 

\MimQX)\^ 3 Q,. \MimQX)\^ ’ 


( 2 . 21 ) 


where the subscript [lc\ indicates that only the leading-colour primitive amplitude (pro¬ 
portional to A^c) is kept. As mentioned above, the DIS-like kinematics associated with this 
antenna function is 

P2 + q^Pi+p-A ( 2 . 22 ) 

withp 2 = P 3 = 0, pf = rriQ, q^ = —Q^ < 0. The corresponding phase space parametrization 
has been derived in [28] in the context of the integration of tree-level initial-final massive 
antennae. It is repeated for completeness in appendix C. 

The phase space integration of the second term in eq.(2.21) is trivial, since the ratio 
jA 42 ^^ ((13)q, 2 g)j^/jA 42 (( 13 )Q, 2q)j^ only depends on Q^, and it can be therefore pulled out 
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of the phase space integral. Thus, only A^{lQ,‘ig,2q) must be integrated over the phase 
space, and this integral is known. 

We must therefore focus on the hrst term in eq.(2.21), for which we follow the same 
procedure as the one described above for ■A^qQgg- The following denominators must be 
considered: 


Di=pI- rr?Q 

T>2 =pI 

D3 = f- m?Q 
D^ = {l-pif 
D5 = (I-Pi -Psf 
D 6 = {I-PI+P2- Psf 

Dj = {p,-qf, (2.23) 


where Di and D 2 are the cut propagators, to Dq are the loop propagators, and Dj = 
—S 23 . We hnd the master integrals depicted in fig.2. 

In addition to these new master integrals, the ultraviolet renormalised integrated an¬ 
tenna -^qQg also contains the known inclusive phase space integral, denoted by /[g] which 
was computed in [24] . This integral is brought about by the phase space integration of the 
mass and strong coupling renormalisation counter-terms. 

The evaluation of the loop integrals Ijg], ^[ 3 , 5 ], d[ 3 , 6 ] aad is rather straightforward. 
It can be done to all orders in e by integrating the underlying loop integrals over the phase 
space. The remaining three master integrals, namely d[ 3 , 4 , 6 ]) -^[ 3 , 4 , 5 ,6] -^[ 3 , 4 ,6,7]; cannot be 

computed in this way. Instead, we calculated them using differential equations. Unlike the 
case of the master integrals needed for ■A^qQgg, the integration constants cannot be in this 
case hxed by imposing regularity conditions on the integrals. This impossibility is due to 
the fact that each integral contains several different powers of (1 — x). Therefore, in order 
to determine the integration constants, we employed in this case independently computed 
soft limits as boundary conditions. For /[ 3 ^ 4 ^ 5 _ 6 ] we used the soft limit of the underlying 
one-loop box given in [42], whereas for d [3 4 g] and /[ 3 , 4 , 6 , 7 ] we calculated the soft limit of 
the one-loop triangle to all orders in e using a Mellin-Barnes representation. More details 
on these derivations are given in appendix C. 

Like the integrated four-parton antenna presented above, -^^Qg has its deepest pole at 
0{\/e^). It can also be written in terms of HPLs and GPLs of trascendentality four and 
three respectively. The HPLs have arguments Xi or xg, and the GPLs have argument xg 
and 1/xj in the weights. The complete expression for -A^q^Qg is enclosed in the ancillary file 
included in the arXiv submission of this paper. Its deepest poles are given by 




(1,Q9 

X 


*7 




(1 -|- Xi) 


- 5(1 - Xi) 


- + + i^o{xi) 


1 


29 -|- 8xgX^ -|- 26xgX^ — IGxgx^ — 52xoXj -|- 5xj 3(1 -|- xj) 


12(1 - XgXi)2 


+ 


4(1 - Xi) 


G(0;xi 
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(d) I[i,e\ (e) -^[3,4,6] (f) -^[3,4,5,61 


(g) -^[3,4,6,7] 

Figure 2: Master integrals for A^^^Qg- Integrals are labelled according to the denominators 
involved (see eq.(2.23)). Bold (thin) lines are massive (massless). The double line in the 
external states represents the off-shell momentum q with = —Q^. The cut propagators 
are the ones intersected by the dashed lines. 


Integral 

Needed in the hard region to 

Deepest pole 

^[0] 

^[ 3 ] 
-^[ 3 , 5 ] 
-^[ 3 , 6 ] 
^[ 4 , 6 ] 
-^[ 3 , 4 , 6 ] 
-^[ 3 , 4 , 5 , 6 ] 
-^[ 3 , 4 , 6 , 7 ] 

e3 

e2 

el 

e2 

e2 

eO 

eO 

eO 

eO 

e-i 

e-i 

e-i 

e-i 

e-2 

e-3 

e-3 


Table 2: Deepest poles and highest order in e needed for each master integral in -A^Qg- 
+(1 + Xi)G{l;xi) - ^ _ ^^G(l;xo) + 2^1 _ 

+(5(1 “ ^ ^ “ ^<^(1; xo) - ^G(l, 1; xq)^ + Vo{xi) + G(l; xq)^ 
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(2.24) 


-2Vi{xi) 


+ 0{e-^) 


3 Top quark pair production at NNLO 

Following the decomposition of the NNLO cross section for top pair production in two, three 
and four-parton hnal state contributions as shown in eqs.(1.2) and (1.4), in this section we 
present the two-parton final state for the coefficient A in eq.(1.5). In particular, we shall 
construct the virtual-virtual subtraction term dd^ nnlo n'^ ^ combination of integrated 

real-real and real-virtual subtraction terms, and mass factorisation counter-terms. We 
will show that these ingredients can be arranged in such a way that dd^ nnlo n'^ 
be expressed in terms of so-called integrated dipoles, which facilitate the cancellation of 
explicit infrared poles in the two-parton final state contribution. 

For massless jet observables computed at NNLO with antenna subtraction, the general 
formalism allowing to write the real-virtual and virtual-virtual counter-terms using inte¬ 
grated dipoles was presented in [18]. There it was shown that massless integrated dipoles, 
it) 

denoted as J 2 with I = 1,2, emerge naturally when rearranging the mass-factorisation 
counter-terms and integrated double real and real-virtual subtraction terms into three dif¬ 
ferent pieces that are free of collinear initial state singularities. It was furthermore observed 
that the colour-ordered J 2 '^’s are related to Catani’s one and two-loop infrared singularity 
operators [43]. 

As we will see below, the virtual-virtual subtraction term dd^- jyz heavy quark 
pair production is similar in structure to the massless case. For the purpose of making this 
similarity manifest, we dehne a new type of massive integrated dipoles J 2 ■ In particular, 
the dipole that we find in the present calculation of this leading colour contribution 
is related to the colour-ordered infrared singularity operator Ig^, defined in [17, 24] and 
recalled in eq.(2.14). 

3.1 NNLO contributions to top quark pair production in the qq channel 

In order to set up our conventions for normalisation factors and matrix elements, we start 
by recalling the leading order partonic cross section for heavy quark pair production. It 
can be written as 

<^^qq,LO=M^^ J d^2{P3,P4-,Pl,P2)\Ml{3Q,AQ,^qAqtJ2^\p3,P4), (3.1) 

( 2 ,^ 

where, d^ 2 {P 3 ,P 4 ',Pi,P 2 ) is the 2 —)■ 2 partonic phase space, and J 2 {P 3 ,P 4 ) is a so-called 
measurement function, which ensures that a pair of final state massive quarks of momenta 
P 3 and Pi are observed. is a colour-ordered and coupling-stripped tree-level ampli¬ 

tude. It is related to the full amplitude through the (trivial) colour decomposition 

~ 3s ~ -^4(3Q)4q, 2g, ig). (3.2) 
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The normalisation factor is 

, ^ ^ (as{n)y Cje)^ (N^ - 1 ) 

LO 2s V 27r J C(e)2 4iV2 

where s is the energy squared in the hadronic center-of-mass frame. Included in A/^lq are 
the flux factor, as well as the sums and averages over colour and spin. The constant C'(e) 
has been defined in eq.(2.3), while C'(e) is given by 


C{e) = 


(3.4) 


providing the useful relation 


= Alias 


C{e) 
271J C(e)' 


(3.5) 


The leading-colour double real corrections to heavy quark pair production are due to 
the tree-level partonic process qq —)• QQgg. They read 


dd 


RR _ 

qq,NNLO,N^ ~ 


^NNLO d$4(p3,P4,P5,Pe; Pi, 1^ (P3,P4,P5,Pe), 


r(4)/ 




(3.6) 


where d<I >4 is the 2 —)• 4 phase space, and \■^^-^^q^^Q^Q^g^gJ' 
coupling-stripped tree-level amplitude normalised to {N^ — !)• 
the overall normalisation A/'nj.^lO’ ''^hich is given by 


is the square of the full 
This factor is included in 


KfRR,qq _ Kfqq 
■'''NNLO ~ ■'^LO 


f as(g) 

V 271 


C(ey 


(3.7) 


The subscript on the matrix element squared indicates that only the terms proportional 
to this colour factor are to be kept. 

Given the fact that the measurement function allows the final state gluons with 
momentum p 5 andpg to become unresolved, i.e. soft or collinear, eq.(3.6) contains infrared 
divergences. These divergences are implicit, in the sense that they only become explicit as 
poles in e after the integration over the phase space is performed. The antenna subtraction 
term that regulates these infrared divergences has been derived in [17]. 

The mixed real-virtual contributions are given by the phase space integral of the in¬ 
terfered one-loop and tree-level amplitudes for the 2 —)• 3 process qq —)■ QQg. They read 


dd: 


RV 

(jg,NNLO,Af| 


= -^NNLO^c / ——d^>3(P3,P4,P5;a;iPl,X2P2)5(l - Xi)(5(l - X 2 ) 


Xi X2 


x2Re (mI , 


M 


Of 


giq2^Q3Qig5 qiq2^Q3Q4,g5 


•/f^(P3,P4,P5), 




(3.8) 


with the normalisation factor given by 

= vfi ( 2 ^)' c(f) = c(t). 


(3.9) 
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The real-virtual contributions in eq.(3.8) contain explicit ultraviolet and infrared di¬ 
vergences as well as implicit infrared ones. The explicit poles in e originate from the loop 
integration in whereas the implicit singularities are due to the phase space 

integration over regions where the matrix elements diverge: the soft limit ps —)• 0 and the 
collinear limits pi||p 5 , P 2 \\P 5 - While the ultraviolet poles are cancelled upon renormali¬ 
sation, we employ a subtraction term dd ^2 to deal with the infrared ones. This 
subtraction term has the twofold purpose of canceling the explicit poles, whose structure 
is well known [44], while simultaneously regularising the phase space integrand in the soft 
and collinear limits. Its explicit construction was derived in [17]. 

Finally, the double virtual contributions have not been derived before and are the 
object of study in the next section. They involve the interference of a two-loop 2 —)• 2 
matrix element with its tree-level counterpart, and a one-loop amplitude squared. They 
can be written as. 


dd 


vv 

(jg,NNLO,iV| 


— ■'''NNLO 


dxi dx2 

Xi X2 


d^2{P3,P4; Xipi, X2P2)S{1 - Xi)(5(l - X 2 ) 






A^\p3,Pi) (3-10) 

JV2 


where the normalisation factor is given by 


■'''nnlo ~ ■'''lo 


a. 


■Xp) 


2tt 


C{er=KNitC{e) 


_ KrRV,qq 
- 


‘'NNLO 


C(e) 


(3.11) 


After ultraviolet renormalisation, these double virtual contributions contain explicit 

infrared poles up to order four, which originate from the loop integration. There are no 

( 2 ) 

implicit poles, as the measurement function J 2 ' does not allow any final state particle to 
be unresolved. As we shall see below, all explicit infrared poles in eq.(3.10) are captured 
and cancelled by those in the subtraction term dd^- ^2 • 

For the two-loop matrix element in eq.(3.10) we employ the analytic results of [33]. 
The square of the one-loop amplitude was calculated in [45]. However, we independently 
computed our own expression, written in terms of GPLs. For the numerical evaluation of 
the HPLs in the double virtual amplitudes we use Chaplin [46], and for the genuine GPLs, 
the GiNaC implementation of [47]. We derived an expansion of both terms in n'^ 

at the production threshold (s ~ 4mg) in order to circumvent the numerical instabilities 
that occur when evaluating the GPLs in this limit. 

Before proceeding to the construction of the double virtual counter-term, it is worth 
commenting on the purpose of the xi and X 2 dependence in eqs.(3.8) and (3.10). This seem¬ 
ingly trivial dependance introduced in dd^j.^j,^j^Q ^2 and dd^^j,^Lo n'^ through delta func¬ 
tions, facilitates the combination with the counter- terms dd^ nnlo ^^^q nnlo n'^ ’ 

which contain integrated antennae and splitting kernels that are non-trivial functions of 
the momentum fractions xi and X 2 - In general, three regions can be distinguished: the 
soft (xi = X 2 = 1 ), the collinear (xi = 1 , X 2 / 1 and xi 7 ^ 1 x 2 = 1 ) and the hard 
(xi 7 ^ 1 X 2 7 ^ 1). The two delta functions in the above equations imply that the real-virtual 
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and the virtual-virtual corrections only contribute in the soft region. Their respective 
counter-terms, on the other hand, contribute in all regions. 


3.2 Structure of the virtual-virtual subtraction term dd^-.TMTr^ at 2 

In general, the virtual-virtual subtraction term denoted as d(Tjjj.^LQ is built from mass 
factorisation counter-terms, as well as integrated double real and real-virtual subtraction 
terms. In the present case, it is given by 


da 


u 

gq,NNLO,Af2 



j ''S,b4 

‘^'^(jq,NNLO,Af2 


/ ‘^'^(jg,NNLO,Af2 ‘^^qg,NNLO,Ar2 
J 2 


f ,-VS,a _ f ,-VS,b _ f ,-.YS,d 

‘^'^<jq,NNLO,Af2 ‘^'^(jg,NNLO,Af2 ^^qq,NNLO,N^ ' 


(3.12) 


The unintegrated forms of all these terms have been given in [17]. /2 ^2 

and contain the new massive initial-hnal integrated antennae presented in 

section 2, namely the tree-level four-parton antenna and the one-loop antenna -^g’qg- 

In general, an additional integrated real-virtual subtraction term, related to a double real 
subtraction term capturing almost-colour-unconnected unresolved limits, and denoted as 
dd^^LQ, should be present in eq.(3.12). Its absence in the present case is due to the fact 
that at the real-real level the partonic process contains only two gluons. 

Following [18] we decompose our virtual-virtual counter-term as 


dd 


u 

qq,NNLO,N^ 


j''U,a 

'^'^qg,NNLO,Af2 


-|- dd 


U,b 

qq,NNLO,N^ 


-|- dd 


U,c 

qq,NNLO,Af2’ 


(3.13) 


where, by construction, each term will be free of explicit initial state collinear poles. In 
order to achieve this, some rearrangements in the integrated real-virtual subtraction terms 
fi d VNNLO,iv2 are necessary. 

We decompose the real-virtual subtraction term j ^2 in two parts: 


f j'-yS,a _ f ,^VS,a,1 , f ,^VS,a ,2 

'a^g5,NNLO,iV2 ~ ‘^^5g,NNLO,Ar2 + ^^qq,Nl<ILO,N2- 


(3.14) 


In the first term, which will be included in ^2 j we group the contributions of the 

form Tg whereas in the second term, which will be a part of dd^'j.^j.^j^Q ^ 2 ; we put 

the subtraction terms of the form . 

The integrated form of the ultraviolet-scale-compensating subtraction term given in 
eq.(8.22) of [17] reads 


i..VS,d 

'a'^qg,NNLO,Af2 


= Ad 


■vv,gg ^2 f da^i dx2 


NNLO 


X1 X2 


d$2(P3,P4;XlPl,X2P2) 








■324) ®2) ®l) 


X jA44(3Q, 4 q, 2g, lg)j J 2 \p3iP4:), 
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(3.15) 


where bo is the Nc coefficient of the QCD beta function at one loop, i.e. 6o = 11/6. The 
square brackets in the second line may be expanded to produce two terms: The terms with 
the pre-factor will be part of ^2 j while those terms originated from 

the factor — 1 in the square brackets with be part of ^2 • 

In order to render each of the three pieces of the virtual-virtual subtraction term in 
eq.(3.13) free of initial state collinear singularities, we must split the mass factorisation 
counter-term ^2 into three parts. This construction is detailed below. 

3.3 The mass factorisation counter-term dd^^i?,. ^ 

g'q',NNLO,A'^ 

In general, for a given partonic process initiated by two partons i and j with momenta 
Pi and p 2 ; the mass factorisation counter-term dd^^^’^^o de included at NNLO at the 
virtual-virtual level reads [18, 22, 23], 


,..MF,2 _ 

*^%',NNLO ~ 


- (£) 


a 




dxi dx2 

Xi X2 


dxi dx2 

Xi X2 




X] ®2)ddfcz,Lo(a;iPi, a;2P2) 


k,l 


dd^i^NLO + dd|]f + / dd® 


kl,NLO 


k,l 


{xiPi,X2P2), 

(3.16) 


where dd^j.,fLQ, dd^^^Q and /y dd|j are the NLO virtual contributions, mass fac¬ 
torisation counter-terms and integrated subtraction terms respectively, and the splitting 
kernels and are given by 

= bki6{l - xi)r[d(x2) -h Sij6{l - X2)r^d(a;^) (3.17) 

rg,L(xi,x2) = Sk^s{l - xi)Tif{x2) + Si, 6 {i - X2)riy (xi) + rid(xi)r[d(x2). (3.18) 

To simplify the construction of the double virtual subtraction term dd^^ can 

be written as suggested in [18] as 

^%ii^UX2) = T[]Ii{xi,X2) - y rjd^;(xi, X2) + ^ 0 (xi, X2), (3.19) 

such that, 

rg.i;(xi,X2) = T^ii\xi)6p6{l - X2)+Tf{x2)6ikS{l - xi). ( 3 . 20 ) 

_^2) 

The reduced kernel PL (z) is related to the usual Altarelli-Parisi splitting functions [48] as 

rlf (^) = (3-21) 

Following [18] the double virtual mass factorisation can be split into three terms, 

1 ^MF,2 _ 1 -^MF,2,a . i '']V[F,2,6 , i ^]V[F,2,c /q oo'\ 

^%',NNLO ~ *^%',NNLO + ‘^%',NNLO + '^%',NNLO’ [O.ZZ) 
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where the individual contributions are given by: 


1 -MF,2,a 
‘^%',NNLO 


a.c 




1 -.MF,2,6 _ 

‘^%',NNLO ~ + 


CX-s 


C(e) 


dxidx2 (1) / V 
dxi dx2 


(^) C'(e)-^d(TM,Lo), (3.23) 


Xl X2 


.(1) 

ij]kl 


-u / fi^s 

'l'^fc«,NLO NLO 


_ /C^ 
V27r 


)<?(') 5 [ryh®r«Jdai,,Lo}. 


,.MF,2,c _ x 2 /■ dxi dj:2 ^(2) 

d%-,NNLO - [211) J Xl X 2 do-fcz^LO- 


(3.24) 

(3.25) 


We apply the generic definitions given above to obtain the mass factorisation counter¬ 
terms dd^ nnlo 7 V 2 '^) leading-colour NNLO corrections to top pair hadro- 

production in the quark-antiquark channel. For ^2 we find 


^^MF,2,a 


gg,NNLO,Af| 


= 

" 27ry 


)cie) 


dxi dx2 


Xl X2 


^qq-,qqi^i^ ^ 2 ) 

bo 


X ( dd^gNLO.TVy - C'(e)yddgg,LO ) (xipi,X 2 P 2 ), 


where dd^ j,^lq is the leading-colour NLO virtual cross section, and 

^il!qq(xi,^2) = 5(1 - X2)rW(xi) -h5(l - Xl)rW(x2), 

with rqg^(x) given in eq.(2.15). 

The second part of the mass factorisation counter-term, dd^ nnlo ’ ^®^ds 


(3.26) 


(3.27) 


, ..MF,2,6 
^^qq,NNLO,N2 


- (^) C{e)]^ [rWgg ® rdy (xi,X2)ddgg-Lo|- 


NLO,TV; 


+ / dd^. 


ijg,NLO,Af, 


(3.28) 


Using the expressions for dd^yo.Ar, and dd®g NLO,iVc given in [17], we get 


dd; 


MF,2,f) 

gq,NNLO 


- [rSgg ® '^ 9 ,Q 5 (u 'Sis)] (a^i, 3^2) - [rjy 0 S24)] (a: 2 , a:i)|dd,q,Lo(a:iPi, 3:2^2)- 

(3.29) 

The convolution of two functions f{xi,X 2 ) and 5 ( 2 / 1 , 2 / 2 ) is defined as 

[/®5](^i,22) = j dxidx2dyidy2fixi,X2)giyi,y2)bizi - xiyi)6{z2 - X2y2), (3.30) 

in such a way that the convolutions appearing in eq.(3.29) take the following forms 
Kq]qq ® ] (xi, X 2 ) = 5(1 - X 2 ) [T^ 0 T^] (xi) + 5(1 - Xl) [T^ 0 T^] (X 2 ) 


- 18 - 








+ 2rW(xi)r«(x2) 


(3.31) 


® sij)] {xi,X2) = <5(1 - X 2 ) [r^ ® ^g,Qg(e, Slj)] (xi) 

+ T[\\x2)AlQg{e, Slj,Xl) 


(3.32) 


[rS„®.45^Qj(e.»23)]te.a'i) =^(1 -a'i)MJ'®.4j,Qj(e.»2j)] fe) 

S%j^X2) 


with 


[f0g]{x)= / dzidz 2 fizi)g{z 2 )S{x - ziZ 2 ). 


(3.33) 


(3.34) 


The convolution of the form [Tgg^ (g) ■^qQg(e, s)] has not been computed before. It 
needed up to finite order and reads 

(8) A^ggie, s)] (x) = s"" X 14 ^6{l - x) + ^^ 0 ( 2 :) 


IS 




^8 

^ - ^5(1 - x) + ^Po(x)) G(l; xo) - ^ ) G(0; x) 

+(1 + x)G(l;x) H- ^ + f I + C 2 ^ <5(1 - x) 


-^'^oix) - 2T>i{x) 


+ 0(0 


(3.35) 


The complete expression can be found in the ancillary Mathematica file attached to the 
arXiv submission of this paper. 

Finally, the mass factorisation counter-term denoted as dd^^^^Q ^2 reads 


dd^^^yp, ^2 = —^c ( ^ ) 0(e 

5g,NNLO,Ar2 ^ \2'K) ^ 


2 1 /=<^,^2 f dxi dX2 =f^(2) 


Xl X2 

T^(2) 


r™;gg(a;i, a:2)dd,q,Lo(a;iPi, 3 : 2 ^ 2 ), (3.36) 


with the leading-colour two-loop kernel Tgg.gg{xi,X 2 ) defined as in eq.(3.20) and 

rfg (^) = ^r(+(x) 

1 r /67 1 2 13 w ^ 1 W ^ N 1 , 2 / N 

--vr -I -ln(x)-ln(x) ln(l — x) -|- - In (x) ]paa(x) 

e[V 72 24 48 ^ ^ 4 ^ ^ ^ ^8 ^ 7 ^ 

+ ^ + -^t) <^(1 - x) 


144 


As usual, pqq is given by 


Pqqix) = 


1 — X 


— 1 — X. 


(3.37) 

(3.38) 
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3.4 The subtraction term ^ 

qg,NNLO,N^ 

The virtual-virtual subtraction term ^2 is given by 


d(T 


U,a 

gq,NNLO,Af2 


- / 


gq,NNLO,Af2 


- / 


— dft 


MF,2,a 


Jl ^ Jl 

\r^V,qq at 2 f *^^2 

■^''NNLO ^''c / 


qq,NNLO,Af2 gg,NNLO,iV| 
d^2iP3,Pi;XiPi,X2P2) 


Xi X2 

'^°,Q3(e>Si3,xi,X2) -rW(xi)(5(l -X 2 ) 
+ ( '^24, X2, Xi) - rW(x2)5(l - Xi) 

A,[lc] 


x<||A^4’""'(3Q,4Q,25,Ig)r--|A^^(3Q,4Q,25,I,)r|>jf (P3,P4) 


(3.39) 


where A 44 ^^‘^^( 3 q, 4Q,2q, Iq) is a shorthand notation for the interference of the one-loop 
leading-colour primitive amplitude A 44 (3q, 4q, 2g, Iq) and its tree-level counterpart. The 

content of the square bracket is manifestly free of initial state collinear divergences as is 
the whole virtual-virtual subtraction term dd^-^T^xT ati■ 

In analogy to the massless case [18], and using a similar notation, we can reformulate 
eq.(3.39) in terms of integrated dipoles. Indeed, we can define the integrated massive NLO 
dipoles as: 


J2(1>'^q) = ^l,Qg{(^^Slj^^Xl,X 2 ) -T^q^{xi)5{l - X 2 ) (3.40) 

J2(2,^q) =^°Qg(e,S 2 j^,X 2 ,xi) -rW(x2)(5(l -Xi). (3.41) 

Using the expression in eq.(2.13) for the pole part of -^qQg one can immediately see that 
these massive integrated dipoles of the form J 2 are directly related to the massive colour- 
ordered infrared singularity operator Iq^. 

Furthermore, the integrated subtraction term ^2 may be rewritten as 

ddJ^’,NNLO,iV2 = -^So^c [ ^^d$2(l>3,P4;xiPi,a:2P2)(j2(i,3Q) + J^(2,4 q)) 

c j X\ X2 X / 

x||Ml’''"’(3Q,4Q,2,-,i,)|2-^|M0(3Q,4Q,2,-,i,)|2| (p3,p4). 

(3.42) 

In general, for an n-jet hadron collider observable computed at NNLO, ddJ^J^Lo contains 
integrated dipoles of three types: final-final, initial-final and initial initial. In our case, for 
the leading colour contributions, the initial-final dipoles defined in eqs.(3.40) and (3.41) 
are sufficient. 
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3.5 The subtraction term dir^-tT.TT ^ 

qg,NNLO,N^ 

For the virtual-virtual subtraction term ^r 2 we obtain: 

qq,NNLO,N^ 

j-U,6 __ _ /*i^VS,6,l _ ,.MF,2,6 

*^^qq,NNLO,N^ ~ / *^^qq,NNLO,N^ / *^^qq,NNLO,N^ ^^qq,NNLO,N^ 

J 2, t/ 1 

= --^NNLO [ ——d$2(P3,P4; a^lPl, 0 : 2 ^ 2 ) { [^^C^q-qq ® (^1> ^ 2 ) 

^ 24 )] (^2? ^ 1 ) \j^qq]qq ®'^q^Qgi^: 

\.'^q,Q9^^^ '* 13 ) ® “^g,Qg('^) '®24)] ) 3^2) 1 1 (^Q, 4q, 2 g, lq)\‘^ J 2 \p^^Pi)- 

(3.43) 


Using eqs.(2.16) and (2.17) it can be noticed that the convolution of two antennae 

S 24 )](xi,X 2 ) is in fact only a product. Furthermore, it can be 
seen that ^2 is not free of explicit initial state collinear poles. In order to rem¬ 

edy this, we define 


1 \],b _ 1 ^ TJ,6 e 1 U 

‘^^g5,NNLO,7V2 ~ ‘^'^qg,NNLO,Af| “ j,jj^LO,iV2 


(3.44) 


with 


-^dd^q nnLO,JV| =-^NNLO [ ^^d^2{P3,Pr,XiPi,X2P2) 


Xi X2 


^ 2 ! •^°,Q9<^^’'®i3) ® ^°,Qg(e,Sl3) {xi,X 2 ) + .4° (c, S 24 ) 0 Qg (e, S 24 ) (X 2 ,xi)| 

X |A44(3q, 4q, 2g, Iq)| J 2 ^{P3iPa)- 

(3.45) 


After this modification, dd^-’^j.^j^j^Q ^2 can be written as a convolution of integrated massive 
dipoles, with no infrared collinear singularities left: 


1 ^ U ,6 _ A rVV,gg a 7'2 

‘^'^gg,NNLO,Ar2 ~ “■'^NNLO 


1 

X - 

2 


dxi dx2 

X1 X2 

rl 


d4’2(4»3,P4;XlPl,X2P2) 


{ [j^(I, 3 q) + J^(2, 4q)] ® [jla, 3 q) + Ji(2, 4 q)] } |M2(3q, 4q, 2,-, i,)|2 ^2) 


X 


2 {P 3 ,P 4 )- 


(3.46) 


It is worth noting that, in general, the terms that we added to dd^’^j^j,jj^Q ^2 in order 
to construct dd^^j,^j^j^Q ^2 are provided by the integrated subtraction term of the form 
dd^^Lo- When those terms are present, which is not the case for our computation as 
discussed previously, dd^),jLo naturally written as a convolution of a sum of integrated 
dipoles without any further additions required. 

Finally, we remark that for the sum of dd^’^j.^j.^j^Q jy 2 and dd^'^j.^j^Q ^2 to remain un¬ 
changed, dd^’'^j.^j^Q ^2 will be redefined as we will see below. 
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3.6 The subtraction term 

qq,NNLO,N^ 

For the integrated virtual-virtual subtraction term ^2 we have; 

J-U.c __ f ,^S,b,4 _ f ,^YS, a, 2 _ f ,^YS,d,2 _ ,-.MF,2,c 

'^^gq,NNLO,Ar2 ~ '^^qg,NNLO,Ar2 ^^qq,N'!<S'LO,N^ ^^qq,N'!<SLO,N^ ^^qq,J<INhO,N^ 

= --^NNLO / — — d$2 (P3, P4; a; iPl, X 2 P 2 ) 

J Xi X2 

^135 a^l; X2^ '^q^Qgg^^'t ^245 ^'2; a;i) 

+dlg,Qg(e, Si3, xi, X 2 ) -F dig Qg(e, §24, 2 : 2 , a:i) 

~l~~ ^ ^2 ^ •^q^Qgi.^J ^lZ^X\iX2') ~\~ ^ ^2 ^ '^24) ^'2; a^l) 

qq,qq{xiTX2^\M.\{2>Q,4:Q,2q,lq)\ \p3,P4)- 

(3.47) 

The expression in the above equation is not free of initial state collinear poles. It 
becomes so, however, once we add back the contributions subtracted from dd^-tT,,. ^ ..2 in 
order to render it free from initial state collinear singularities. Indeed we find that 

‘^'^gg’,NNLO,Ar2 ~ '^'^q^,NNLO,Af2 + ^^^qq,N'NhO,N2 (3.48) 

has a similar structure than the one found for ddJ^J^Lo ™ massless case, and allows for 
the definition of new initial-final massive integrated dipoles: 

J2(i) 3 q) = dig Qgg(e, Si3, ^1, ^ 2 ) + dl^ Qg(e, Si3, xi, X 2 ) -I- — 

-rgg(xi)5(l-X2)-^ dl°Qg(e,Si3)(8)^°Qg(e,Si3) (xi,X2) (3.49) 

J2(2; 4q) = dig Qgg(e, ^24, X2, Xl) -|- ^g Qg(e, S 24 , X2, Xl) -|- — 

-rjg(a:2)5(l-a;i)-^ .A°Qg(e,S24)«>dl°Qg(e,S24) (a;2,a:i). (3.50) 

In terms of these massive integrated dipoles we have 

ddJJ,NNLO,7V2 = --^NNLO^c / ^ ^d4>2(p3, P4; , X 2 P 2 ) 

X (ji(I, 3 q) + Ji(2, 4 q)) \MI{3q, 4q, 2,-, lq)f (^ 3 ,^ 4 ). (3.51) 

As in the massless case, one can see that the integrated dipoles present here are related 
to integrated subtraction terms which involve genuine NNLO objects. This feature can be 
regarded as an additional check on the correctness of the construction of the subtraction 
terms at real-real, real-virtual and virtual-virtual levels in this NNLO computation. 
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3.7 Explicit pole cancellation 


With the explicit expressions included above one can show analytically that each of the 
virtual-virtual subtraction terms dd^-^T.TT ^ ..o i dd^-and is free of 

gg,NNLO,A'^’ qq,[\NLO,N^ qq,NNLO,N^ 

explicit initial state collinear poles and that the sum of those building blocks is such that 

Voles (^d(T^^j^LO,Af2 “ dd^ j<fNLo,Ar 2 ) “ (3.52) 

This is an ultimate check that demonstrates that we have correctly implemented the sub¬ 
tractions terms at real-real, real-virtual and virtual-virtual levels for the leading-colour 
contributions to qq —)■ QQ -|- X at NNLO. 


4 Heavy quark contributions 

In this section we present the heavy quark contributions to the top pair production cross 
section in the quark-antiquark channel computed at O(a^). They are proportional to the 
colour factors NhNc (leading-colour) and Nh/Nc (subleading-colour), and correspond to 
the coefficients and Gh in eq.(1.5). 

Both Ffi and receive contributions in the two and three-parton final states, and 
require, in both final states, subtraction terms to cancel explicit and implicit infrared diver¬ 
gences. For all loop matrix elements involved we perform the ultraviolet renormalisation in 
the so-called decoupling scheme, with the gluon wave functions and the heavy quark mass 
and wave functions renormalised on-shell, and the strong coupling ag renormalised in the 
MS scheme with Ni active flavours. The two-loop amplitudes employed in our calculation 
[33], were originally computed with the strong coupling renormalized in the MS scheme 
with Np = Nfi + Ni active flavours, and therefore required a conversion to the decoupling 
scheme. 

In this section we will explicitly construct the two and three-parton contributions to 
the NfiNc and N^/Nc colour factors, and we will show how the amplitudes computed 
in [33] in the full-flavour scheme can be converted to the decoupling scheme via a finite 
renormalization of Ug. 


4.1 Real-virtual contributions 

In the antenna subtraction framework, the three-parton heavy-quark contribution to top 
pair production at NNLO in the qq channel is given by 



d(T 


RV 


d(T 


T 

gg,NNLO,Vh 


(4.1) 


Following the colour decomposition of the real-virtual matrix-element as presented in [29] 
can be written as 


^^qq,NNLO,Nh 


= -^NNLO [ ^'^d<^3{P3,Pi,P5;XiPi,X2P2)S{l - Xi)(5(l - X 2 ) 


Xi X2 


xiNr 


|M^’‘^J(3q, 5g, i,;; 2,-, 4 q)| 2 + i,;; 2,-, 5„ 4 q)P 
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(4.2) 


1 


5„ 4q; ; 2,-, i,)\^ + (3q, 4^;; 2,-, 5„ i,)p 


-2|A^^’^^<(3Q,4Q,25,i„5^)P 


«/2^^(P3,P4,4'5), 


where we have used the following shorthand notation: 

| Af ^’["'(...)|2 = 2 Re ( M ^’'"'(---)- M 5 (---)*) 


(4.3) 


with and At 5 ’^^^ being respectively the tree-level and the Nh coefficient of the one-loop 
amplitude associated to the process qq —)• tig. In addition, in eq.(4.2) we have defined the 
one-loop amplitude AI 5 ’^^^( 3 q, 4q, 2q, Ig, 5..y) as 


(3q, 4q, 2,-, 5^) = (3q, 5„ 1,;; 2,-, 4q) + (3q, 1,;; 2,-, 5„ 4q) 

= -^ 5 '^ ^(3q, 5g, 4 q; ; 2g, ig) -|- ^ (3 q, 4q; ; 2g, 5^, ig), (4.4) 




where the gluon is photon-like. 

After UV renormalisation in the decoupling scheme dd^f.,fj,^LQ is free of explicit 
e-poles. It contains, however, infrared implicit poles: The phase space integral is divergent 
due to integration over the soft limit P 5 —)• 0 and the collinear limits pi | |p 5 and ^ 2114 * 5 ) where 
the matrix elements are singular. The subtraction term dd)^ in eq.(4.1) regularises 

these implicit divergences of the phase space integral. We shall present it below. 


4.2 Real virtual subtraction term dd^ 

Both the leading-colour (NhNc) and subleading-colour (Nh/Nc) real-virtual subtraction 
terms only contain terms commonly known in the antenna subtraction literature as dd^®’“^^^. 
They are built as products of tree-level antennae and reduced one-loop matrix elements 
squared. Unlike in the most general case, there are no integrated double real subtraction 
terms, no mass factorisation counter terms dd'^^’^^ ^or subtraction terms involving one- 
loop antennae. The absence of one-loop antennae in ddA 

the Nh coefficient of the antennae Ag, which would be in principle be required, vanishes 
in the decoupling scheme. It is also worth mentioning that in the decoupling scheme the 
reduced one-loop matrix element needed in dd, i.e. the part of the one-loop 
amplitude for the process qq —)• ti, is finite. 

Explicitly, we hnd. 




= -^NNLO f —— d4>3(p3,P4,P5; XiPi, X2P2)S{1 - Xi)(5(l - X 2 ) 


xlNr. 


Xl X2 

, 1 , 


A0(3Q,5„I,)|Ml’''^‘((35)Q,4Q,2,-,i,)|2jf (P3g,P4) 


+A0(4Q,5„25)|Ml’'"'(3Q,(45)Q,2,-,i,)|Vf (p3,P45) 


1 


2A0(3Q,53,2,-)|M^’^'^'((35)Q,4Q,2,-,I,)|24"^(p3g,P4) 
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(4.5) 


-2A0(3Q,5„ig)|A^4''^’((35)Q,4Q,2,-,l,)|Vf (p35,P4) 
—2A3(4q, 5g, 2g)|A44’^ ^(3 q, (45)q, 2q, lq)| J 2 \p3,P^) 
-^°(3Q,5,,4Q)|All’['^]((3~5)Q,(4~5)Q,2,-,i,)|Vf 
-^°(2,-,5„i,)|A^^’''^'(3Q,4Q,2,-,l,)|2jf (P3,P4) 


Since, as we mentioned above, |A^4’^^^(3 q, 4q, 2^, lq)P has no e-poles, the subtraction term 
in eq.(4.5) does not contain any explicit infrared singularities. 

Being a pure subtraction term, dd^^ nnlo Nh be added back in integrated form 

to the cross section at the two-parton level. The required A-type integrated massive and 
massless antennae are all known [24, 27]. 

4.3 Double-virtual contributions dd^^js^LQ 

In this section we present the heavy quark NNLO two-parton contributions to qq —?■ tt in 
the decoupling scheme, which are given by 

/ [do’^^NLO.Af^ “ NNLO,Af^] • (4-6) 

4d<I>2 

We shall focus in particular on the construction of the double-virtual counterterm 
b'5'^NNLO Vh’ iiito a leading-colour (iV/iW) and a subleading-colour (Nh/Nc) part. 
The contributions from the double-virtual matrix elements can be written as 

dd^g^NNLO,iVh = ^NNLO^h [ ——d4>2(p3, ■ ■ • ,P4; Xipi, X 2 P 2 ) 

J X\ X2 

x<5(l - xi),5(l - X2)|M1(3 q,4q, 2,-, 4"\p3,P4) , (4.7) 

with A/’nj^]('o given in eq.(3.11), and with the abbreviation 


(4.8) 


Nh 


For the part of the two-loop matrix element in eq.(4.8) we employ the analytic results 
of [32] converted to the decoupling scheme via a finite renormalisation of ag that will 
be described below. The “one-loop squared” term has been computed analytically in 
[45]. We re-derived it ourselves, also analytically, and use our own result in our event 
generator. We further compared these results with those provided by Roberto Bonciani in 
both renormalisation schemes and found full agreement. 

After UV-renormalisation, dd^^j^LQ contains explicit infrared poles, that are can¬ 
celled by the virtual-virtual subtraction term dd^ nnlo Nh derived below. 







4.4 Virtual-virtual subtraction term 

Following the decomposition of the mass factorisation counter term and dd^ ex¬ 

plained in section 3, we find that dd^j,^jvjLQ only contains terms denoted there as dd^’“, 
involving only the mass factorisation counterterm dd^^^^^Q and the integrated form 
of the real-virtual counter term derived in eq.(4.5). 

The mass factorisation counterterm is given by 


T^MF, 2 ,a _ \ f 

a^gg,NNLO,Ar;, - y — 


Nh 


Nr_ 


dxi dx2 
Xi X2 


X2)dd^g^NLO,ivJa;ipi, X2P2) 


(4.9) 

with Fgg-qg given in eq.(3.27), and where dd^j^jLQ^^ the Nh coefficient of thes the NLO 
virtual cross section. 

We decompose the double virtual counter term into a leading and a subleading-colour 
part: 


dd 


u 

qq,NNLO,Nh 


dd 


u 

qq,NNLO,NhNc 


-|- dd 


U 

qq,T<INhO,Nf,/N,- 


(4.10) 


The leading-colour part is given by 


dd^,NNLO,AfhAfc 


= -A/'nnLO ^hNc / — —d4>2(p3,P4; XiPi,X2P2) 
J Xi X2 


si 3 > a:i, X 2 ) - rW(xi)(5(l - X 2 ) 


+ ('®24> X2,xi) - rW(x2)(5(l - xi) 


m 




(4.11) 


The terms present in the square bracket are manifestly free of initial-state collinear diver¬ 
gences as is the whole double virtual subtraction term dd^ iV^Afc' Using the same 
integrated massive NLO dipoles as in section 3.4 for the leading-colour double virtual 
counter-term contributions, we can rewrite eq.(4.11) as: 

dd^5,NNLO,iV;,JV, = —d4>2(p3,P4; , X 2 P 2 ) 

J X\ X2 

X (ji(I, 3 q) + Ji(2, 4 q))(3q, 4q, 25 , ig)| (p 3 ,P 4 ). (4.12) 

The subleading-colour heavy quark counter term ddU n^/n '^^S'ds, 

.rVV,qqNh f da^l dX2 , . , 

d%-,NNL0.7V;,/7V, = -^NNLO]^ J —Pi] Xipi, X2P2) 

X ^2.Aq Qg(e, S23, X2, Xl) -|- 2jAq Qg(^€-, Sj4, Xl, X2) ‘^•^q,Qg{^J '^13’ ^1’ ^2) 

2.Ag Q^(e, 524, ^2, Xl^ '^34, Xl^ X2) '^qq^gi^^i '^12? ^1? ^2) 

+rW(xi),5(l - X2) + r«(x2)5(l - Xl ))%, i,)|2jf (p3,P4). 
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(4.13) 


This subtraction term can also be written in terms of integrated dipoles. In addition to 
the initial-final massive dipoles dehned in eqs.(3.40-3.41), this reqnires a new final-final 
massive dipole defined as 

3l{lQ,KQ)=Al^Q{e,siK), (4.14) 

as well as a massless initial-initial dipole dehned in [18] and given by 

3l{tq, %) = A^gg^g{xi,X2) ” rW(xi)5(l - X 2 ) ” rW(x2)5(l - Xi). (4.15) 

In terms of these three type of integrated dipoles, the virtual-virtual subtraction term 
d<7S,NNLO,iV;,/7V, takes the following form: 

J^U _ ArVV,(jq f dxi dx 2 . 

^^qq,NNLO,Nh/Nc, — ~-^NNLO at / —d4>2 (p3, P4; , X 2 P 2 ) 

-‘''C J XI X2 

X (2ji(2, 3 q) + 2ji(I, 4q) - 23l{l, 3q) - 23^(2, 4q) - J^(3q, 4q) - 3^(1, 2)) 

x|Ml’''^'( 3 Q, 4 Q, 25 ,i,)| 2 jf (p 3 ,P 4 ). ( 4 . 16 ) 

The pole parts of each of the integrated dipoles in eqs.(4.12) and (4.16) is given by 
a specihc colour-ordered infrared singularity operator given in [17]. In terms of those 
operators we hnd that the pole part of the entire double-virtual heavy-quark counter term 
reads: 


'Poles i d(T^q NNLO,Af,,Afc 


= -NhM, 


^^^q,NNLO,Nh/Nc 


W,qq f dXl dX2 


NNLO 


Xi X2 


d^2{P3,PPXiPi,X2P2)S{l - Xi)(5(l - X 2 ) 


8 iVjg(e, ^13) - ^ ( 4 lg(e, S13) - 41 ^( 6 , S23) + 4 ? ^12) 

Re{Rh)\Ml’^''\3QAQ, %, 4)P 4"\p3,P4). 




(4.17) 


As mentioned above, the renormalized one-loop amplitude is hnite. It is fur¬ 

thermore proportional to its tree-level counterpart: 


-^4 ^ ^ (3q, 4q, 2g,lg) — Rh A44(3q, 4q, 2g,lg), 


(4.18) 


with the factor Rh in the equation above given by 

Rh = -l{8-3^^)-^(3-^^)H{0-,x)-in^(3-(3‘^)+O{e) (4.19) 

and with (3 = — 4771 ^/ 512 . Using this expression for the one-loop amplitude, we can 

also express the pole structure of the “one-loop squared” contributions in dd^^j^LO 
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given in eq.(4.7) in terms of infrared singularity operators. For the leading-colour part we 
find 

4lg)(e, sis)Re{RH)\Ml{3QAQ, 2g, (4-20) 

and for the subleading colour part, we get, 

Voles (3q, 4q, 2g-, (3 q, 4q, 2,-, = 

2(4lg)(e,si3)-4lg)(e, 523)+4^(6, si2) + lJ5^(e,S34))Re(R,.)|A4°(3Q,4,3,2, 

(4.21) 


Taking the pole part of the integrated subtraction terms given above (where we have used 
momentum conservation to set 524 = si 3 and 514 = 523 ) enabled us to verify the pole 
structure of the two-loop amplitudes in the decoupling scheme. 

We conclude this section by showing how the one and two-loop amplitudes computed 
in [33] in the full-flavour scheme can be converted to the decoupling scheme with Ni active 
flavours. This conversion can be achieved with the well known finite renormalisation of as 
(see e.g. [49, 50]) 






27r 


27r 


2 

6 + 0(«s) , 


(4.22) 


with denoting the strong coupling in the full theory, and the strong coupling 

in the decoupling regime. The renormalisation constants are given by 


si ~ 2 


6 = 
with 


1|+£^(5C^ + 3Cf) 




4 —+ 3Cp)L^^^ — Nfi ^ 

(47r 


2 \ e 


r(l + e) -e 


m 


Q 


o-«7b 


(4.23) 

»■» 


(4.25) 


We can use eq.(4.22) to relate an amplitude Ai computed in the full theory to its 
counterpart in the decoupling scheme, denoted as Ai. We have, 

2 ^ 


. , (W)\ 

Ai = (dvraf-)) WfW + ^ + 



a. 


27T 

m' 


271 


+ 0{a 


277 






+ 0(a: 


?) , 


(4.26) 
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such that for the one and two-loop amplitudes we find respectively, 

;^(l) (4_27) 

(2) ^ ^ + (^ - 2 (6)^) (4.28) 

5 Numerical Results 

Together with the double real and real-virtual leading-colour contributions derived in [17], 
the two-parton channel presented in section 3 completes the calculation of the coefficient 
A in eq.(1.5). We implemented this coefficient along with Fi, F^, Gi, Gh, Hi, Hh and Hih, 
in a Monte Carlo parton-level event generator based on the set up of eq.(1.2). While the 
inclusion of the coefficients Hx is rather straightforward, as they are all infrared finite and 
only receive contributions from double virtual matrix elements, the coefficients Fx and 
Gx are non-trivial. Fi and Gi were presented in [29], and TX and Gh were derived in this 
publication in section 4. 

Using an adaptation of the phase space generator employed in the NNLO di-jet calcu¬ 
lation of [51], our program, written in Fortran, provides full kinematical information on 
an event-by-event basis, thus allowing for the evaluation of differential distributions for top 
pair production including exact NNLO results for the quark-antiquark channel. Naturally, 
the calculation of total cross sections is also possible. 

In this section, we present numerical results for LHC and Tevatron. For LHC, with 
^/s = 8 TeV, we show differential distributions in the transverse momentum of the top 
quark pip, the top quark rapidity y^, the invariant mass of the tt system and its 
rapidity y**. For Tevatron we show differential cross sections in pp, and m^i, as well as 
in the absolute value of the top quark rapidity jy*j. In a separate sub-section we present 
our results for the forward-backward asymmetry Aps, differential in the rapidity difference 
= \y^ ~ y^\ well as in and pp. 

All the results presented below include all partonic channels at LO and NLO. At NNLO, 
the contributions to the quark-antiquark channel from the colour factors N^, NiNc, Ni/Nc, 
NhNc, Nh/Nc, Nf, NlNiNh and NhN^ and Nh/N, are included. We use the pole mass 
of the top quark m* = 173.3 GeV, and the PDF sets MSTW2008 68cl. The factorisation 
and renormalisation scales are set equal to the top quark mass pp = pp = "mt throughout. 
Where provided, scale variations correspond to the range mt/2 < rrit < 2mt. 

The total hadronic cross section in the quark-antiquark channel is presented in table 3 
for Tevatron, and LHC with center-of-mass energies 7, 8 and 14 TeV. From the comparison 
with the full NNLO calculation of [52] it can be seen that, due to the omission of subleading- 
colour pieces at O(a^), our results consistently overestimate the NNLO correction to the 
total cross section, indicating that the combined contribution of the coefficients G and 
E in eq.(1.5) is negative and non-negligible. It should be noted that the size of these 
subleading-colour terms decreases with the hadronic center-of-mass energy, as can be seen 
from the fact that the discrepancies between or results and those of [52] at NNLO are 5% 
for Tevatron, 3.25% for LHC 7, 3% for LHC 8 and 2.25% for LHC 14. 
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LO 

NLO 

NNLO (LC) 

NNLO (Ref [52]) 

Tevatron 

6.069 

5.845 

6.425 

6.118 

LHC 7 

29.76 

28.21 

29.51 

28.59 

LHC 8 

37.60 

35.67 

37.31 

36.22 

LHC 14 

91.22 

85.91 

89.96 

87.99 


Table 3: Contributions to the total hadronic cross section from the quark-antiquark chan¬ 
nel at LO, NLO and NNLO at different collider energies. Cross sections are in [pb], with 
fJ-R = ^^F = Frit- Our NNLO predictions are contained in the third column, labelled as 
NNLO (LC), as they are dominated with the leading colour (LC) part. 

5.1 Differential distributions for LHC 

In fig. 3 we present differential distributions for top pair production in pp collisions with 
y/s = 8 TeV at LO, NLO and NNLO, together with the corresponding /c-factors. As 
expected, the impact of the NNLO corrections in the qq channel on LHC cross sections is 
in general mild, given the dominance of the gluon-gluon initiated process in this scenario. 

From the ratios NNLO/NLO in the lower panels of figs, (a) and (c), it can be seen 
that the NNLO corrections decrease the ptp and distributions over the entire spectra 
considered. The decreases range between 3% and 7%, being in both cases more pronounced 
in the tails of the distributions. 

In the distributions in and y**, given in figs, (b) and (d), it can be seen that in 
comparison with the full NLO result, NNLO corrections in the quark-antiquark channel 
shift the cross sections downwards in the central region by 5%. In the very forward and 
backwards ends of the spectra the impact of these corrections is more substantial, causing 
an upwards shift of 12% in the y* distribution, and over 50% in the y** distribution. 

All plots in fig. 3 show a slight reduction in the scale uncertainty when the NNLO 
corrections to the qq channel are included. 

In order to assess the relative size of the contributions from the different colour factors 
included at NNLO in our calculation, in fig. 4 we show the breakdown of the NNLO cor¬ 
rections to the qq channel as functions of pip and y*. We find that the leading-colour piece, 
proportional to N^, contributes most significantly, followed by NiNc, which was calculated 
in [17]. The contributions from all other colour factors are very small. 

To further disentangle the light and the heavy quark contributions, in fig 5 we show 
only the fermionic contributions, omitting the colour factor. 

5.2 Differential distributions for Tevatron 

In fig. 6 we present differential distributions for top pair production in pp collisions with 
^/s = 1.96 TeV at LO, NLO and NNLO, along with the corresponding ratios. Given 
the dominance of the quark-antiquark channel in this scenario, the impact of the NNLO 
corrections included in our calculation is more important here than in the LHC distributions 
presented above. In all cases we find good agreement with experimental data, and a 
significant reduction in the scale uncertainty at NNLO. 
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= l73.3.Gey 

/^i? ~ /^F ~ '^hop 

MSTW2008 68clpdf 



(a) 


(b) 



Figure 3: Differential distributions for LHC {^/s = 8 TeV) in: (a) top quark transverse 
momentum pip, (b) top quark rapidity y^, (c) invariant mass of the tt system (d) 
rapidity of the tt system y**. The NNLO contributions included are in the quark-antiquark 
channel only. Renormalisation and factorisation scales are set equal pR = pp = and 
varied as mt/2 < p < 2mt. 


As can be seen from the distribution in pip given in fig.(a), in the lower part of the 
spectrum, NNLO corrections introduce a 10% shift with respect to the NLO prediction. 
This shift decreases as pip increases, becoming negligible in the tail of the distribution. In 
the distribution in |y*| given in fig.(b), NNLO corrections to the qq channel introduce an 
increasing shift of up to 30% with respect to the NLO result. A significant reduction in the 
scale uncertainty is also observed. In fig.(c) we show the inclusive cross section for top pair 
hadro-production as a function of the invariant mass of the top-antitop system. NNLO 
corrections, in this case, cause a positive shift over the entire spectrum ranging from 15% 
near the production threshold to approximately 5% in the tail of the distribution. 

It is well known that QCD corrections to the y* distribution in pp collisions cause a 
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(b) 


Figure 4: Contributions to the NNLO QCD corrections to ppiqq) tt + X from the 
different colour factors included in our computation. 




(a) (b) 


Figure 5: Contributions to the NNLO QCD corrections to pp{qq) tt + X from fermionic 

colour factors. 


shift towards the forward region. This fact can be observed in our result shown in fig.(d), 
where it can also be seen that NNLO corrections modify this shift in comparison to the 
NLO result. As we will see in the next section, this will have an effect in the NNLO result 
for the forward backward asymmetry. 

In fig. 7 we show the breakdown of the NNLO corrections to the y* and distributions 
into colour factors. The same pattern found in fig. 4 can be observed, with the colour 
factor contributing most significantly, followed by NiNc- As before, in fig. 8, we omit the 
leading color factor in order to further examine the size of the different fermionic colour 
factors. 
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Figure 6: Differential distributions for Tevatron in: (a) top quark transverse momentum 
pip, (b) absolute value of the top quark rapidity \y^\, (c) top quark rapidity y*, (d) invariant 
mass of the tt system rritt- The NNLO contributions included are in the quark-antiquark 
channel only. Renormalisation and factorisation scales are set equal pn = pp = and 
varied as mtopl‘2: < y < Experimental data points from CDF and DO are taken 

from the results in the i -t- jets channel presented in [3] and [4] respectively. 


5.3 The forward-backward asymmetry at Tevatron 


The forward backward asymmetry in hadronic top quark pair production is an observable 
that measures the difference between forwardly and backwardly produced top quarks. It 
is defined as 


a(Ay^) - cT(AyL^) 
(T(Ayf^) -h o-(Ay“) 


(5.1) 


with Ay** = y* — y* and Ay^ = y(±Ay**). 

At leading-order App vanishes given the forward-backward symmetry of the Born 
matrix elements. At higher orders, however, loops and real emissions in the quark-antiquark 
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Figure 7: Contributions to the NNLO QCD corrections to ppiqq) tt + X from the 
different colour factors included in our computation. 




(a) (b) 


Figure 8: Contributions to the NNLO QCD corrections to pp{qq) tt + X from fermionic 

colour factors. 


and quark-gluon channels cause the (anti)top to be more likely produced in the hemisphere 
defined by the direction of the incoming (anti)quark. Due to this fact, which was first 
observed in [53], a non-vanishing Aps m pp collisions is predicted starting at next-to- 
leading order. 

Much attention was dedicated to the forward-backward asymmetry after measurements 
by the CDF collaboration showed a pronounced discrepancy with the Standard Model NLO 
prediction [5, 54]. A recent publication by the DO collaboration [55] showed measurements 
that differ from those by the CDF collaboration, and are in agreement with all SM pre¬ 
dictions. In [14], the NNLO QCD corrections to App were calculated for the first time, 
yielding a significant reduction in the scale uncertainty and confirming the agreement be- 
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tween the SM prediction and the latest result from DO. 

Here we present our results for the forward-backward asymmetry at NNLO, computed 
in the so-called unexpanded form: 


^FB 


a^AaNLO + oigAaNNLO + C>(«s) 
ajcTLo + a^cTNLO + ajaNNLO + C>{a^) 


(5.2) 


In fig. 9 we show Afb as a function of = \y^ — y^\ as well as in and p^. In 

the first two cases we find agreement with the DO measurements, as well as a substantial 
reduction in the renormalisation and factorisation scale dependence at NNLO. The change 
from NLO to NNLO in the scale variation bands of the distribution is rather unusual, 
with a reduction in the scale dependence in the first bin, and an increase in all others. 

We compared the results in hg.9 for the differential asymmetries with [14]. Although 
at NNLO our computation only includes the quark-antiquark channel at NNLO (and it 
does not include all colour factors), when this channel is dominant, the two results are in 
agreement. Our results appear to differ at the edges of the distributions. In those regions 
the top quark pair are produced mostly asymmetrically where we expect the quark-gluon 
channel, not included in our computation (but included in [14]), to play an important role. 


6 Conclusions 

We presented the computation of the virtual-virtual two-parton final state contributions 
to top pair production in the quark-antiquark channel proportional to the leading colour 
factor N^. Together with the three and four-parton contributions derived in a previous 
publication [17], this enabled us to complete the NNLO corrections to top pair hadro- 
production in the quark-antiquark channel for the phenomenologically most important 
colour factor. 

We derived the subtraction term needed at the two-parton level using the antenna 
subtraction method extended to deal with massive coloured final states. This required 
the integration of new massive tree-level four-parton and three-parton one-loop antennae, 
which we presented here for the first time. We showed that the virtual-virtual counter-term 
has a similar structure to that found in the massless case [18], and can be expressed in terms 
of integrated massive dipoles whose pole part is related to the massive infrared singularity 
operators of [44] . We furthermore demonstrated that the explicit poles of the double virtual 
matrix elements can be analytically cancelled against those of the corresponding counter¬ 
term, providing a crucial check on the correctness of our calculation. In addition, we 
presented the computation of the heavy quark contributions related to the colour factors 
NhNc and Nh/Nc- In particular, we presented our results using analytic expressions for 
the loop amplitudes in the so-called decoupling scheme, which were not available as such 
in the literature. 

Together with the NNLO corrections proportional to the number of light quark flavours 
Ni obtained in [29], as well as other easily added fermionic colour factors, we implemented 
the leading-colour and heavy quark NNLO contributions to qq ^ tt + X completed in this 
paper, in a fully differential Monte Carlo event generator. This program allowed us to 
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Figure 9: Differential Aps in: (a) absolute rapidity difference |Ay**| = \y^ — y*|, (b) 
invariant mass of the tt system and transverse momentm of the tt system The 
NNLO contributions included are in the quark-antiquark channel only. Renormalisation 
and factorisation scales are set equal yp = yp = y and varied as mtopl"^ < Ai < “^rntop- 
Experimental data points from CDF and DO are taken from the results in the i + jets 
channel presented in [54] and [55] respectively. 


produce several differential distributions for top pair production with NNLO corrections 
in the quark-antiquark channel for Tevatron and LHC energies. In particular, we obtained 
distributions in the top quark transverse momentum pip and rapidity y*, as well as in the in¬ 
variant mass and rapidity of the top-antitop system and As expected, we found that 
at the LHC, the NNLO corrections to the qq initiated process are not phenomenologically 
significant, except in the very forward and backward regions of the rapidity spectrum. For 
Tevatron, on the other hand, we found that these corrections are important and drastically 
reduce the scale dependence in all distributions considered. 

In addition, we evaluated the NNLO corrections to the differential forward-backward 
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asymmetries at Tevatron and found good agreement with the results from DO [55]. In the 
regions where the contributions from the quark-gluon channel are expected to be small, we 
also found agreement with [14]. 
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A Master integrals for A^qgg 

In this appendix we collect the six master integrals found in the IBP reduction of the 
integrated antenna A® Qgg . The Laurent expansion of these integrals can expressed in 
terms of HPLs with arguments Xi or xq and GPLs of argument xq and weights involving 
1/xj. We keep the overall factor 

A^xo 4 = C'(e)^ + w-q) ^^(1 — 2:) (A.l) 

unexpanded, in order to facilitate the introduction of the master integrals into the inte¬ 
grated antennae and the subsequent expansion in distributions. 

For the integrals /[g] I[-s] we have the following all order expressions [24]: 

^[0] = Nxo4 exp^^T'® (g2 + 

X 2 E 1 fl-e,2-2e,4-4e; / ~ ) (A.2) 

V 1-XiXoJ 


/[_ 8 ] = Nx04 


rHi-e) 

2r(4 - 4e) 


exp^^T'B ^ 2^2 2+2f. 


-2+2e 


X 2 E 1 1 - e, 2 - 2e, 5 - 4e; 


1 - 

1 - XjXg 


(A.3) 


These all order expressions can be expanded up to the required order in e using the 
Mathematica package HypExp [56]. The expanded forms written in terms of HPLs and 
GPLs are included in the ancillary file attached to the arXiv submission of this paper. 
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The remaining master integrals are computed using differential equation techniques. 
We had already calculated /[ 4 ] in [25]. Written in terms of Xi and xq it reads: 


^[4] = 


^X04 

Xo 


xoG(0; Xi) - (1 - xo)G(l; xq) + (1 - xo)G ( —; xq ) + e 

Xi 


2'k^xq 


+5xoG (0; Xi) - 5(1 - xq)G (1; xq) - (1 - xq)G (0; x*) G (1; xq) + 5(1 - xo)G —; xq 

\Xi 

+3(1 - xq)G (0; Xi) G ( —; Xo ) + SxqG (0,0; x*) - 2G (0,1; xq) + 2G ( 0, —; xq 

\Xi ) \ Xi 

+4xoG (1,0; Xi) + (1 - xo)^ (1,1; xq) - 3(1 - xo)^ ^1, 


Xi 


+(1 - xo)G —, l;xo + (1 - xo)G —, —;xo 


1 1 


Xi Xi 


+ O(e^) 


(A.4) 


The masters /[4 gj, .^[ 4 , 5 , 8 ] j .^[ 4 , 7 , 8 ] new. The first terms in their Laurent expansion 
are given by 

I[4:,8] = Nxoa + mQ)”^ Xi |g( 0, l;xo)G(0;xi) - G{0;xi)G l;xo^ 


vr 


1 


Xi 


-—G —;xo + G(l;xo)G(0,0;xi) - G(0,0;xi)G —;xo - G(l;xo)G(l,0;Xj) 


1 


1 


Xi 


1 


+G(1, 0;xi)G I — ;xo ) — G I 0,0, — ;xo ) — G \ —,0,l;xo ) + G I —,0, — ;xo 


Xi 


1 


Xi 


1 1 


Xi Xi 


2 2 

+—G(l; Xo) + G(0,0,1; xq) + —G(0; Xj) + G(0,0,0; Xi) — G(0,1, 0; Xi) + 2^3 + 0(e) 


(A.5) 


-^[ 4 , 5 , 8 ] — -^X 04 (Q^ + nrg) 


2^-2 Xf 


1 1 
1 — Xi 1 36 ^ 


G(0;xi) + -G(l;xo) + 


1 

+ - 
e 


-G(l; xo)G(0; Xi) + -G(0; Xi)G I —; xo + -G I 1, —; xo + -G I —, 1; xo 


Xi 


1 


1 


1 


1 


Xi 


1 „ / 1 1 


Xi Xi 


1 


57r 


24 


+ xG* ( —, —;xo ) — -G(l,l;xo) + 3G(0,0; Xj) —— 


6 

+ 4-2 


18 


+ 0(e° 


(A.6) 


-^[4,7,8] — -^X04 (Q^ + "mg) Xj <| — ^ — ^G(l; Xo) + - 

1 


1 


Xi 


77r2 

- ^ - G(0,l;xo) 
3d 


+G I 0, —; Xo ) + 2G (l,0;xi) + -G(l,l;xo) — Gl 1, —;xo 


1 


Xi 


+ G(e° 


(A.7) 


B Massive phase space for A^ ggg 

In order to provide a boundary condition for the system of differential equations involving 
the integrals 7[4,8]) %, 5 , 8 ]! ^[ 4 , 7 , 8 ] given in the previous section, we computed the soft 

limit of 7[4,7,8] via direct evaluation. This required a parametrization of the phase space 
associated to the DIS-like kinematics 

P 2 + g ^ Pi+P 3 +P 4 (B.l) 
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with P 2 = pi = p\ = 0, pf = rriq, = —Q^ < 0 . As far as we know, this phase space has 
not been derived previously in the literature. We present it below. 

The starting point is 


d^>3 = 


1 


(27r)2(i 3 

1 


d<ipid"*P3d"*P4(5+(pf - rnQ)5 ^{pI)6^{pI)5^‘^\pi +P 3 +P4 -P2 - g) 


P3 




4(27r 


l2d-3 


E3 E4 


5+((P3 +P4 -P2 - qf - mq) 


(B.2) 


where we have used the momentum-conserving delta function to integrate out pi. Now, we 
parametrise the momenta as 


P 2 = E2{l,0(d-2), 1) 

q = {Ecm - T' 2 , 0 (rf_ 2 ), -E 2 ) 

P 3 = £^ 3 ( 1 , 0 (d- 3 ), sin 01,008 01 ) 

P 4 = £^4 (1,0 ((i_ 4 ) , sin 02 sin 03 , sin 02 cos 03 , cos 02 ) 
Pi=P 2 + q-P3 -P 4 , 


(B.3) 

(B.4) 


such that, 

j d^-^p3 = 
f d^-^p4 = 


dEs £|-2 J d(cos 0 i) sin'^'-^ 0 i J df 2 d -2 

dE 4 E^~^ J d(cos 02 ) sin*^”"^ 02 J d(cos 03 ) sin*^”^ 03 J 414^-3. (B.5) 


The phase space measure takes the form 

^d-2^d-3 


d4>3 = 


4(271)2*^ 3 J 

1 d—3 


X ^ 3^4 

with the solid angle given by 


d£' 3 d£ 4 d(cos 0 i)d(cos 02 )d(cos 63 ) 

1 d — 4: 


sin 01 sin 62 

— 


sin''-5 03. 


nd/2)- 

We can now change variables to the following invariants; 


(B. 6 ) 

(B.7) 


523 = 2p2 • P3 = 2£2£3(1 “ COS0i) 

524 = 2p2 • P3 = 2£2£4(1 “ COS 62 ) 

534 = 2p3 • p4 = 2^3£4(1 — sin 01 sin 62 cos 63 — cos 0 i cos 62 ) 

513 = 2pi • P3 = 2p2 ■p3 + 2q-p3- 2p3 • P4 = 2£3£cm - S 34 

514 = 2pi • P4 = 2p2 • P4 + 2^ • P4 - 2p3 • P4 = 2E4,Ecra “ S 34 . (B.8) 


The corresponding Jacobian is 

9(si3, S14, 523, 524, 534 ) 


J = 


det 


d{E 3 , £ 4 , cos 01 , cos 02 , cos 03 ) 


1 -1 


32£2^£|£|£| sin 0 i sin 02 


(B.9) 
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Note that the energy E 2 is fixed such that 2 p2 • q = and S 12 is not an independent 

integration variable. It is given by 


S 12 = 2pi ■ P 2 = {E1^ + 0^) — S 23 — S24- (B.IO) 

Using the above expressions, the phase space measure in d = 4 — 2e dimensions reads: 
f 2(47r)-^+2^ 


r'i—4H-2e /* 

d4>3 = ~r(i _ 2g) ^ J dsi2dsi3dsi4ds23ds24ds34(-A4)“^/^"" 

x 0 (—A 4 ) 5{eI^ — rriQ — S 13 — S 14 — S 34 ) 5{E^^ + — S 12 — S 23 — 'S 24 ) 


(B.ll) 


where the dependence on S 12 has been introduced through a delta function, and the Gram 
determinant is given by 


A 4 = A(si2S34, 'S13'S24, 'S 14 S 23 ) + 4mQS23S24'S34 


with 


A(x, y, z) = x'^ + y‘^ + — 2xy — 2xz — 2yz. 

We now define the following dimensionless variables: 


m, 


and 


1 


U 2 = 


Z23 = 


u{E‘^m - ^Q. 

1 

- m^Q, 


uo = 


S 12 


Q 


Ecm - ’ 


u = 


EL + Q^ 

E^m - 


(B.12) 

(B.13) 

(B.14) 


1 


1 


U3 = 


Sl3 


•523 


{E'^m - m^q) 

1 

2^24 - -71Z3 -9T S24 


Zi4 = 


{E^m - mq'. 
1 


Sl4 


u{E‘^m - "ig) 


234 = 


(E'^m - m^q. 


S 34 


with 


(Ecm - rnq) = — —(Q^ + m^) 

Xi 


where 


Xi = 


<3^ + "ig _ <3^ + ^q 


and 


2P2W {E^^ + ml) (3^+"ig' 

In terms of the dimensionless variables Zij the phase space reads. 




(B.15) 

(B.16) 

(B.17) 


2(47r)-4+2^ , 
d$3 = ^ ^1^2 


{Ecm + Q 


2\l-2e 


d2;i2 d2;i3 d2:i4 d2;23 d2;24 d2:34 


r(l-2e) 

X (—A 4 ) 0(—A 4 ) 5(1 — 2:13 — zi 4 — Z 34 ) 5(1 — 2:12 — 2:23 — 224 ) 


with 


A 4 = \{zi 2 Z 34 , 2:13^24, 2142:23) + 4rio223224234- 


(B.18) 

(B.19) 
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We use the delta functions to integrate out zi 2 and Z 34 , and reparametrise the remaining 
integration variables as 


2^13 = XI 

X2(l - Xi) 


Z24: = 


Z23 = 


(1 + Uo) 

XiX3(l + '«o -X 2 (l-Xi)) 


(1 + uo)(no + Xi) 

Zl4 = Z ^4 + X4(^14 - Z 44 ), 


(B.20) 


where are the roots of A 4 viewed as a function of 2 : 14 . The Jacobian for this reparametri- 


sation is 


, , , rxi(l-Xi)(l+ «o - (1-Xi)X2)l 

= I - (H-„o)Vo + Xl) - I 


dxidX2dx3, 


(B.21) 


and the integration region gets mapped to the hypercube 0 < Xi < 1 - We finally arrive at 
r- 4 + 2 e 


d$3 = 


TT 


i 28 r(i - 




a- 2 e^-l+ 2 e 


/ 


l-2e 


xi(i - Xl) 

(1 - Xo) + Xl (1 - Xi) 


dxidX2dx3dx4 

11 / 2 -e 


X2(1-X2)J [X3(1-X2)J [X4(l-X4) 

1 -1+e 


(B.22) 


In the soft limit x* —)• 1, after the trivial integration over X 4 is performed, our phase 
space parametrisation completely factorises as, 


dd) 


3,soft 




1 


+ rriQ)' '^"{1-Xi)^ ^'(l-xo) 


\-2+2e 


X / dxidx2dx3 X 


Xi(l - Xl) 


l-2e r 


X2 (l-X2 ) X3(l-X3) 


(B.23) 


and the invariants take the following simple form 

■S 12 —^ + ihq) 

si 3 ^ + m|) (1 - Xi)xi 

■S23 —^ + in'Q)xiX 3 

•524 —^ + "iQ)X2(l — Xl) 

514 {Q"^ + mg)- -^214 


Xi 

2 I ^2 \ (1 “ Xi) 


Z 34 


534 {Q + iHq) 

Xi 

(1 — x) r 

5234 -(<3^ + "^q) ^\ _ [xiX3 + (1 - Xl)X2 


(B.24) 


Using eqs.(B.23) and (B.24) we can find explicit all order expression for all master 
integrals that do not involve the invariants S 14 or S 34 . As a boundary condition for the 
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integrals found in the reduction of we only needed to compute the soft limit of 

-^[ 4 , 7 , 8 ]) which reads 


Asoft) 

[4,7,8] 


3-^2 (1, —2e, 1 — e; 1 — 3e, 1 — 2e; 1) 


(B.25) 


C Master integrals for 

In this appendix we collect the seven master integrals found in the IBP reduction of the 
integrated antenna -A^q^Qg- The Laurent expansion of these integrals can expressed in terms 
of HPLs with arguments Xi or xq and GPLs of argument xq and weights involving I /xj. 
The phase space associated to the DIS-like kinematics of this antenna, namely 


P2 + q^ Pi+P2, 


(C.l) 


with P 2 = pi = 0, PI = rriQ, = —Q^ < 0, was derived in [24], It is given by 


(dvr) 


— 2 +e 


J Y(^i _ g) "xl{l-Xi) "(1-XiXo) 

with 


-l+e 


"^0=712 




Q + m' 


Xi = 


Q 


Q'^ + ttIq "'■ 2 p2-q 


dyy "(1-y) S 
(C.2) 


(C.3) 


The inclusive phase space integral, denoted as /[q], is obtained by performing the y 
integration in eq.(C.2): 


-^[0] = iQ'^ + mg) " (dvr) 


^-l+e 


r(l-e) 


X-(l-Xi)^ ^^(1-XjXo) 


-l+e 


2r(2 - 2e) 

With this parametrization of the phase space, the invariants can be expressed as 


(C.4) 


Sl 2 = {Q +"7g) 


2 f 1-y- XjXo + Xiy 


023 -«3 


(1 - Xi)y 


x(l - XjXo) 

The third invariant S 13 can be obtained using momentum conservation. It reads 

/ 1 — Xj 

S13 = S 12 + S23 -Q’^ -mQ = {Q'^ + mg)' 


(C.5) 


(C. 6 ) 


The evaluation of the master integrals Ip], ^[ 3 , 5 ], L[ 3 , 6 ] ™d /[4 6 ] is straightforward. 
It can be done to all orders in e by employing an all order representation of the tadpole 
and bubble integrals involved, and integrating over the phase space in eq.(C.2). The loop 
integrals that we need are 


Tad(m|) = J 


d'^l 


1 


(27r)'^ P — m? 


- jc r(i - 2e)r(-l + e) , 2 

r 2 {i - e)r(i + 3) (“O' 


(C.7) 
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d'^/ 1 


= l- 


Cr 


Bub(/c2;0,0) J ^27rYP {I-k)^ ^e{l-2e) 


i-e) 


, 2 \-e 


(C.8) 


Cr r(l- 2 e) , 2 ,-. ^ f 

= 1— --r (m^) 2-^1 ( e, 1 , 2 - e;— 


with 


e(l-e)r 2 (l-e) 

1 r2(l-e)r(l + e) 


m 


(C.9) 


( 47 r) 2 -^ r(l - 2 e) 

Integrating these loop integrals over the antenna phase space in eq.(C.2) we find 


(C.IO) 


/[ 3 ] = j d<^Xij, Tad(m^) 

= iVxi 3 S, (e^ + ml) 


e(i_e)(i_ 2 e) 


(C.ll) 


I[ 3 , 5 ] = j Bub(si3 + 0 , m^) 

Nx13 S, 


e(l - 2 e )2 
„l+e 


°”(2")r(i-^)r(i + 2^) ^3.(j _ ^^,2-4., 1 _ ^.^^,-2+2. 


2 r(l + e) 


+Xj +''(1 - Xi)^ ^^{ 1 -xqY ^{l-XiXo) 2'+'= 2^1 1^1,6, 2e; Y 

%,6] = [ d$x,.^,Bub(-Q^0,m2) 


1 — X 


XXo 


(C.12) 


Nx 13 Se 


x-(l-Xi)^ ^"(1-XiXo) (1 - e,e,2 - e;xo) (C.13) 


eil-e){l-2e) 

%, 6 ] = J Bub(-S23, 0 , 0 ) 

with the overall factors Nxi 3 and given by 

Nxi3 = iC{efiQ^ + ml)-^^ 

and 


S, = Cr, 


o<^lE 


(C.14) 


(C.15) 


(C.16) 


2 r(l-e)- 

The remaining three master integrals, namely -I[ 3 , 4 , 6 ]) -^[ 3 , 4 , 5 ,6] -I[ 3 , 4 , 6 , 7 ]) cannot be 

computed in this way. This would require all order expressions for the triangle and box 
integrals depicted in fig. 10 and given by 

/ dYl 1 

Tr^772-2T7i-1271-1-12 

(277)"* [P - m^){l -piY[l-pi +P2 -P3) 
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Figure 10: Triangle and box one-loop integrals 


Box(m^, 0, 0, S 13 + -S 23 ; 0,0, 0) = 


dH _ 1 _ 

{2-kY {P - rrpQ){l - piY{l - Pi - pzY{l -P 1 +P 2 - PzY '' 


(C.18) 


which are not known. 

We computed these mixed one-loop phase space integrals with differential equations, 
using their soft limits, independently evaluated to all orders in e, as boundary conditions. 
For /[ 3 ^ 4 ^ 5 ^ 6 ] ''^6 used the soft limit of the underlying one-loop box given in [42], whereas 
for /[ 3 , 4 , 6 ] s-iid /[ 3 , 4 , 6 , 7 ] we calculated the soft limit of the one-loop triangle to all orders in 
e using a Mellin-Barnes representation. 

We express the integrals in terms of the regular functions as in eq.(2.9) in order to 
facilitate the expansion of in plns distributions. The decomposition of .^[ 3 , 4 , 6 ]) -^[ 3 , 4 , 5 ,6] 

in) 

and /[ 3 , 4 , 6 , 7 ] into regular functions Ra ’ reads 

/[3,4,6] = iVxi3((l - ^)“'^^[|,4,6] + (1 - ^)"''^[3,4,6] + (^ " ^)"''^[3,i,6]) (C-19) 


'[ 3 , 4 , 6 , 7 ] 


^X13((l - x) 


-2ep(2) 

[ 3 , 4 , 6 , 7 ] 


+ (1 - x)-^^R 


—3e t?(3) 


[ 3 , 4 , 6 , 7 ] 


+ (1 - 


—4e d(4) 


[ 3 , 4 , 6 , 7 ] 


(C.20) 


^[3,4,5,6] = ^X13(l - x) '^^^[3^4_5 


(C.21) 


(n) 

The first few terms in the Laurent expansion the Ra '’s in the above equation are given by: 


R 


( 2 ) 


1 


[ 3 , 4 , 6 ] 


1 


G(l;xo) -G (—■xq 


1 


Xi 


1 


4G ( 0, —; xo ) - 4G (0,1; xg) 

Xi 


1 


R 


(3) ^ _ 

[ 3 , 4 , 6 ] Q2 _|_ jyi2 


Q 


Xi< - 


Xi 


4^2 


—G (1, 1; Xq) -l- 3G I 1, — ;xo) + G( — ,l;xo) — 3G I —, —; xq 


1 


1 1 


Xi Xi 


+ 0(e°) 


G(l;xo) - G ( —;xo 

Xi 
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